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Abstract 

Using a self-consistent quark model for nuclear matter we investigate variations 
of the masses of the non-strange vector mesons, the hyperons and the nucleon in 
dense nuclear matter (up to four times the normal nuclear density). We find that 
the changes in the hadron masses can be described in terms of the value of the scalar 
mean-field in matter. The model is then used to calculate the density dependence 
of the quark condensate in-medium, which turns out to be well approximated by a 
linear function of the nuclear density. Some relations among the hadron properties 
and the in-medium quark condensate are discussed. 
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I INTRODUCTION 



One of the most interesting future directions in nuclear physics may be to study how nu- 
clear matter properties change as the environment changes. Forthcoming ultra-relativistic 
heavy-ion experiments are expected to give significant information on the strong inter- 
action (i.e., QCD ) in matter, through the detection of changes in hadronic proper- 
ties m, 0, |] . In particular, the variations in hadron masses in the nuclear medium have 
attracted wide interest because such changes could be a signal of the formation of hot 
hadronic and/or quark-gluon matter in the energetic nucleus-nucleus collisions. Several 
authors have recently studied the vector-meson (u, p, 0) masses using the vector dom- 
inance model [|], QCD sum rules |J and the Walecka model (QHD) || [F], |8|], and have 
reported that the mass decreases in the nuclear medium. There is also a proposal to look 
for such mass shifts at CEBAF [[|. 

In the approach based on QCD sum rules, the reduction of the mass is mainly due 
to the four-quark condensates and one of the twist-2 condensates, (q^^D^q). However, it 
has been suggested that there may be considerable, intrinsic uncertainty in the standard 
assumptions underlying the QCD sum-rule analyses [JlO | . On the other hand, in hadronic 



models like QHD ||11|| , the main reason for the reduction in masses is the polarization of 



the Dirac sea |], [7], ||, where the antinucleons in matter play a crucial role. However, 
from the point of view of the quark model, the strong excitation of nucleon-antinucleon 
pairs in medium is difficult to understand. It is clear that these two mechanisms are quite 
different. 



Several years ago Guichon |12| proposed an entirely different model for nuclear mat- 
ter, based on a mean-field description in which quarks (in nucleon bags) interact self- 
consistently with a and u mesons. The model has been refined by Fleck et al. [l3j and 



the present authors. It provides a natural explanation of nuclear saturation and the right 
magnitude of the nuclear compressibility. We have used this model to investigate the 
nuclear structure functions |TJ]], the properties of both nuclear and neutron matter jjTJJ 
and the Okamoto-Nolen-Schiffer anomaly and isospin symmetry breaking in matter |16| . 



We argued that the response of the internal structure of the nucleon to its environment is 
vital to the understanding of, not only physics with momentum transfers of several GeV 
(eg., deep- inelastic scattering), but also physics at scale of a few MeV (eg., a violation of 
charge symmetry in nuclear medium). In Ref.||15||, the relationship between the Guichon 
model (alias the quark-meson coupling (QMC) model in our papers) and QHD has also 
been clarified. Even though this model is extremely simple, the insights gained from it 



may help to reveal the essential physics []I7| . Here we use the Guichon model to inves- 
tigate variations of various hadronic properties, as well as the density dependence of the 
quark condensate in nuclear matter, as functions of the nuclear density up to four times 
normal nuclear density. 

This paper is organized as follows. In Sec. II, the QMC model is introduced and 
applied to calculate the modification of the properties of the nucleon in nuclear matter. 
The masses of the vector mesons (a;, p) and the hyperons (A, E, E) in matter are studied 
in Sec. III. We find that the change in the hadron masses can be described in terms of 
scalar mean-field values in medium. Some relationships among the hadron masses are 
given. Next, in Sec. IV, we use this model to study the density dependence of the quark 
condensate, which may be linked to a diverse range of nuclear phenomena. The density 
dependence of the quark condensate can be well fitted by a linear function of the nuclear 
density. Then, some connections among the hadron properties and the quark condensate 
in medium are discussed. Sec. V contains our conclusions. 



II MATTER PROPERTIES IN DENSE MEDIUM 

II. 1 The quark-meson coupling model 

The QMC model treats nuclear matter as a collection of (nucleon) MIT bags, self- 
consistently bound by the exchange of scalar (a) and vector (u, p) mesons. We assume 
that nuclear matter with N ^ Z is uniformly distributed, and that the mesons can be 
treated in mean-field approximation (MFA). Let the mean- field values for the a, cu (the 
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time component) and p (the time component in the third direction of isospin) fields be 
a, uo and b, respectively. The quarks in a static spherical bag interact with those mean 
fields. The Dirac equation for a quark field, ip q , in a bag is then given by 

[z 7 • d - (m 9 - V a ) - 7 °(K, + l -r z v p )]^ q = 0, (1) 

where V a = g%a, V u = g q u and V p = g^b, with the quark-meson coupling constants, g 9 ,, 
g^ and g q . The bare quark mass is denoted by m q and r 2 is the third Pauli matrix. The 
normalized, ground state for a quark is then given by 



where 



e q = Q q + R(V U ±\V P ), for P ) quark, ( 3 ) 



N- 2 = 2R 3 j 2 (x q )[n q (n q - 1) + Rmy2]/xl (4) 
/?, = ^(fi,-i&n$)/(fi, + i*m*), (5) 



with f2 9 = ^/a;;j + (Rm*) 2 and the quark spinor. The effective quark mass, m*, is 
defined by 

m* = m q - V a , (6) 

for both u and d quarks. The linear boundary condition, jo(x q ) = f3 q ji(x q ), at the bag 
surface determines the eigenvalue, x q . 

Using the SU(6) spin-flavor wave function for the nucleon, the nucleon energy is given 
by E^ ag + 3V W ± ~V P for , where the bag energy is 

EL = £ '"'"°"^ + (7) 



with I? the bag constant and a phenomenological parameter accounting for a multitude 
of corrections, including zero-point motion. Here n q is the number of quarks in the nucleon. 



To correct for spurious cm. motion in the bag [y| the mass of the nucleon at rest is taken 
to be 



MN = J{EZ a )*- y En q {x q /R)*. 



The effective nucleoli mass, M N , in nuclear matter is given by minimizing eq.(|8|) with 
respect to R. 

Table 1: B 1 ^ and z N for some bag radii (mo = 5 MeV). 



Ro(fm) 


0.6 0.8 


1.0 


B 1 ^ (MeV) 


187.7 157.2 


136.1 


z N 


2.038 1.640 


1.169 



To see the sensitivity of our results to the bag radius of the free nucleon, R Q , we choose 
the current quark mass, m q = m (= m u = m^) = 5 MeV, and vary the parameters, B 
and zn, to fit the nucleon mass at 939 MeV with Rq (= 0.6, 0.8, 1.0 fm). The values of 
B 1 ^ and zn are listed in Table []]. 

For infinite nuclear matter we take the Fermi momenta for protons and neutrons to 
be k F . (J = p or n). This is defined by pj = k F J(37i 2 ), where pj is the density of protons 
or neutrons, and the total baryon density, p_e, is then given by p p + p n . 

Since we want to calculate the variations in the masses of the u and p mesons, we 
suppose that both mesons are also described by the MIT bag model in the scalar mean- 
field. To fit the free masses, m w = 783 MeV and m p = 770 MeV, we introduce new 
z-parameters for the two mesons, z u and z p . We then find that z w { p ) = 0.7840(0.8061), 
0.4812(0.5160), 0.1198(0.1680) for R = 0.6, 0.8, 1.0 fm, respectively. Unfortunately, in 
this model it is hard to deal with the density dependence of the a meson in medium because 
it couples strongly to the pseudoscalar (tt) channel, which requires a direct treatment of 
chiral symmetry in medium ||19|| . That is beyond the scope of this study. Though one 
might expect the cx-meson mass in medium to be less than the free one (see, eg., Ref. ||), 
we shall keep the free value, 550 MeV, even in medium. 

The io field is now determined by baryon number conservation as Co = g^pB/m^, 
where g u = 3g^ and is the effective omeson mass, and the p mean-field by the 
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difference in proton and neutron densities, P3 = p p — p n . On the other hand, the scalar 



mean-field is given by a self-consistency condition (SCC) |TT], T^|. Since the p field value 
is b = 5^3/ (2m* 2 ), where g p = g q p) the total energy per nucleon, E tot , can be written 



E, 



tot 



E 



dk\/Mtf + k 2 + 



m 



o- ^2 



2pi 



r + 



2m* 2 8m* 2 p# 



Pi 



(9) 



where m* is the effective p-meson mass in medium. Then, the SCC for the a field is given 
by 



a 



(27r) 3 m2 
Using Eqs.(0) and (||), we find 



E 

j=p,n ■ 



dk 



Mi 



Mf + k 2 



'dM 



N 



da 



(10) 



R 



'dM 



N 



da 



-g a C N {a), 



where g a = 3g% and Cn is the quark-scalar density in the nucleon 

a, 



C N (a) 



^bag 



M 



N 



1 - 



S N + q 



E N 



with 



S 



N 



R 



drip q i> q 



n q /2 + Rm*{Sl q - 1) 



(12) 



(13) 



Q q (Q q - l)+Rm*/2 

The value of the quark-scalar density is about 0.4 at ps = 0, and it decreases monotoni- 
cally to about 0.2 at the normal nuclear density, po = 0.1 7fm~ 3 . The dependence of the 
scalar density on the bag radius is not strong (see Ref.|T5|). 



II. 2 Coupling constants and nuclear matter properties 

We determine the coupling constants, g\ and g 2 , so as to fit the binding energy (—16 MeV) 
at the saturation density, po, for symmetric nuclear matter. Furthermore, the p-meson 
coupling constant is used to reproduce the bulk symmetry energy, 33.2 MeV. The values of 
the coupling constants are listed in Table [|. In the last two columns, the relative changes 
(with respect to the values at zero density) of the bag radius and the lowest eigenvalue are 
shown. The present model gives a good value for the nuclear compressibility - around 200 
MeV. If we take a heavier current-quark mass, the effective nucleon mass and the nuclear 
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Table 2: The coupling constants and calculated properties of equilibrium matter at the 
saturation density (mo = 5 MeV). The effective nucleon mass, M^, and the nuclear 
compressibility, K, are quoted in MeV. 



R (fm) <£/4tt gl/An g 2 p /4-K M* N K 



SR 



Sx 



XQ 



0.6 



18.79 



1.326 4.923 838 223 -0.03 -0.07 



0.8 



20.63 



1.016 5.014 850 205 -0.02 -0.09 



1.0 



21.09 0.871 5.045 855 196 -0.01 -0.12 



Figure 1: Scalar mean-field values. The dotted, solid and dashed curves are, respectively, 
for #0 = 0.6, 0.8 and 1.0 fm. 

compressibility are, respectively, reduced and enhanced. The strength of the scalar mean- 
field, V a , in medium is shown in Fig.|l]. At small density it is well approximated by a 
linear function of the density (see also Eq.(40) in Ref . [|15| ) : 



Figure 2: Effective nucleon mass in symmetric nuclear matter. The curves are labelled as 
in Fig.l. 

The dependence of the effective nucleon mass on the nuclear density is shown in FigfJ. 
It decreases as the density goes up, and it behaves like constant at large density. We find 
that the effective nucleon mass at small density is approximately given by 




(14) 




(15) 
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which is identical to the model independent result derived using QCD sum-rules 



Figure 3: The ratio of the axial-vector coupling constant of the nucleon in medium to 
that in free space. The curves are labelled as in Fig.l. 



Figure 4: The ratio of the magnetic moment of the proton in medium to that in free 
space. The curves are labelled as in Fig.l. 



In Fig.^] the density dependence of the axial- vector coupling constant of the nucleon 
is illustrated. It decreases as the density increases. On the other hand, the magnetic 
moment of the proton in symmetric nuclear matter increases with density as shown in 



Fig.^. The calculation of these quantities is standard [21] (cm. and recoil corrections are 
not included |22|). It is rather easy to understand why g\ and ^ N tend to decrease and 
increase, respectively, with density: in the bag model these quantities involve integrals over 
the upper and lower components of the quark wave function. Because the attractive scalar 
potential effectively decreases the quark mass it makes the solution more relativistic - i.e., 
the lower component of the wave function is enhanced. This fact leads, respectively, to 
the decrease and increase of g\ and fi* N in medium. Given the simplicity of this argument 
we find the conclusion that the magnetic moments of the other hadrons should increase 
quite compelling. 

At small density, we can expand g\ and /j,* n in terms of ps- If we take the quark to 
be massless for simplicity, and ignore the dependence of g\ and p* N on the change of the 
bag radius in medium, they are expressed as 

g\ 4x2 - 12X 2 , + 10x - 3 , , ( p B \ 

9a 2x5^X0 -l) 2 \poJ 
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and 

~ 1 - 4xg ~ 16xg + 17xq - 6 ~ i + .1 , (17) 

/ia 2x (x - l) 2 (4x - 3) CT \Po/' 

where x = 2.04, i? = 0.8 /m and Eq.(|TJ]) is used. These formulae can roughly reproduce 
both the ratios at small density. Furthermore, we have calculated the charge radius of the 
proton in medium, and we find that the change is within a few per cent even at ps ~ 4po- 

Here we must add a caution that the contribution of meson exchange currents (MEC) 
to g% p* N , etc. would be considerable in medium, and, hence, we should treat the whole 
problem including MEC to get the final results for the changes in these quantites. This 
is obviously beyond the scope of the present work. 

Finally, we record that the variations of the above nucleon properties depend only 
very weakly on the proton fraction, f p , which is defined as p p /pb, because the p meson 
does not play a role in determining the nucleon structure in medium (see Eqs.(Q) ~ (§)). 
Further investigations concerning the nucleon properties, the equation of state for matter, 



etc. can be found in Ref. M, ITR 116 
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Ill HADRON MASSES IN MEDIUM 

Now we are in a position to present our main results. In Fig.|5|, the ratio of the effective 



Figure 5: The ratio of the vector- meson mass in symmetric nuclear matter to the free 
mass. The curves are labelled as in Fig.l. 



vector- meson mass, m*, to that in free space is shown as a function of the density. Since 
the difference between the effective u- and p-meson masses at the same density is very 
small, we show only one curve for both mesons in the figure. As the density increases 
the vector- meson mass decreases (as several authors have previously noticed @, H |, 
and seems to become flat like the effective nucleon mass. The mass reduction can be well 
expressed by a linear form at small density: 

<~1-0.09(^V (18) 
m v \ p J 

The reduction factor, 0.09, is somewhat smaller than that predicted by the QCD sum 
rules ||. In this model the reduction in the mass is basically caused by the attractive 
scalar mean- field in medium, and this origin is clearly different from that in QHD || || , 
in which the essential mechanism is the vacuum polarization due to the excitation of 
nucleon- antinucleon pairs in medium. 

It is possible to calculate masses of other hadrons in this model. In particular, there 
is considerable interest in studying the properties of hyperons in medium - eg. A, S and 
S. For the hyperons themselves we again use the MIT bag model, including the cm. 
corrections (see Eq.(H)). We assume that the strange quark does not couple to the a 
meson in MFA, and that the addition of a single hyperon to nuclear matter of density ps 
does not alter the values of the scalar and vector mean-fields. The mass of the strange 
quark, m s , is determined so as to reproduce the mass splitting between the nucleon and 
the average of A and E hyperons, Mn av (= 1154 MeV), in free space. We then find m s 
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= 355.5, 358.1, 354.0 MeV for R = 0.6, 0.8, 1.0 fm, respectively. Using these values of 
the strange-quark mass, we have calculated the average mass of A and E, , and the 
mass of H, MS, in symmetric nuclear matter. 

Figure 6: The ratio of the average mass of A and E in medium to that in free space. The 
curves are labelled as in Fig.l. 



In Fig.[], the average mass, M|j- , in medium is illustrated. The density dependence 
of Mt, is very similar to the cases of the vector meson and the nucleon. The ratio is 
again well described by a linear function at small density: 

~ 1-0.08 ( — ) . (19) 
M Hav \Po) 

Note that the reduction factor in the mass formula is almost the same value as in Eq. (|18"D . 



Figure 7: The ratio of the S mass in medium to that in free space. The curves are labelled 
as in Fig.l. 



The effective S-hyperon mass in medium is shown in Fig.[7|. The ratio again behaves 
like the case of M^ av , but an approximate form for MS/M^ at low density is given by 

— - ~ 1 - 0.04 [ — ), 20 

where the reduction factor is about half of those in Eqs.(|i~8"D and dl9|). 

We summarize the effective masses of the hadrons in medium in Fig.|| As one can 
see from the figure, the reduction in both the vector-meson mass and the average mass of 
A and E are about twice that in the S-hyperon mass over a wide range of ps, while the 
reduction in the nucleon is about three times that in the H. One can also see this fact 



from the reduction factors in Eq. (|T5|) and Eqs.([T8|) ~ (|20|). It is rather easy to understand 
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Figure 8: The ratios of the effective masses of the hadrons in symmetric nuclear matter 
to those in free space (Ro = 0.8 fm). The dashed, solid, dotted and dot-dashed curves 
are for H, (A + S)/2, the vector (u, p) meson and the nucleon, respectively. 



why such a relationship holds among the hadron masses: the nucleon consists of three 
non-strange quarks, which couple to the scalar mean-field in medium, while the vector 
mesons, and the A and £ hyperons, involve just two non-strange quarks. The H hyperon 
consists of only one non-strange quark and two strange quarks. Since the strange quark 
does not feel the scalar mean-field, its mass is not altered in medium. The change in the 
hadron mass is therefore roughly determined by the number of non-strange quarks, no, 
and the strength of the scalar mean-field, V a . We can then find an approximate form of 
the hadron mass in symmetric nuclear matter: 



hadron ~ 1 - const x n V a (MeV), (21) 

^hadron 

where const ~ 3.2 x 10 _4 MeV _1 . Note that this formula is not too sensitive to the bag 
radius, and that it can reproduce the hadron masses reasonably well over the range of 
pB up to about 3po- From Eq. (|2~ID we expect the following relations among the hadron 
masses: 



(Mk) ~ (Mi) 3 (<) ~ (Mi) ~ (Mi) ~ (Mi) 2 (22 ) 

\M N ) ~ \M E ) ' \mj ~ \M A J ~ VM s y ~ \M E J ' 1 ; 

Furthermore, using Eqs.(|l^) and (p^7|), g% p* N and m* at small density could be linked as 

(<w^w^y\ (23) 

\m v J \g A ) \m ) 
Finally, we record that the proton-fraction dependence of the effective masses we have 
studied in this section is again very weak. 
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IV QUARK CONDENSATES IN MEDIUM 



Having shown that the QMC model provides an interesting description of the hadron 
mass in nuclear matter we now apply it to study the behavior of the quark condensates in 
medium. The quark condensates are very important parameters in the QCD sum rules, 
and it is believed that they are linked to a wide range of nuclear phenomena, including 
the effective hadron mass in medium || ^ . 



The difference of the quark condensate in matter, Q(pb), and that in vacuum, Q(0), 
is given through the Hellmann-Feynman theorem fl5], [16], |l7|, [24j| : 

1 dS 



Q( PB )-Q(Q) 



2 dm q 
OS dM 



N 



dM^ dm q 



+ E 



dm* dm„ 

j=mesons J 1 



dS dm* . , 88 da; 

— - + J2— + 



dgj dm q 



(24) 



where £ = p B E tot . 

Using the SCC, one finds 



88 



dM 



N 



dM 



N 



and 



drrin 



dm* 



Cn(o-) V 9u 



3C N (a) 1 - 



dK 
dm,, 



2CJa) 1 



dV a 



dm q \ dm q 



(25) 
(26) 

(27) 



where v — u or p, and C v is the quark-scalar density in the vector meson given by Eq.(|T2"D 
(using the variables for the vector meson instead of those for the nculeon). Then, using 
the parametrization for the derivative of the cx-meson mass with respect to the quark 
mass |2l, ^ = ^Tm^' wnere a ^( = ^rn q C N (0)) is the nucleon a term f25|, and the 



Gell-Mann-Oakes-Renner relation, the ratio of the quark condensate in medium to that 
in vacuum is written by a rather lengthy formula: 



Q(Pb) 

Q(Q) 



l 



m. 



— ) {g a a) + A x a 2 + A 2 p\ + A 3 p 2 3 

9a J 



(2* 
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where m n is the pion mass (138 MeV), f n the pion decay constant (93 MeV) and the 
coefficients, Ai ~ A 3 , are given by 

A^CMO^-U^) 2 ^, (29) 
M N 6 \g a J dm q 

A 2 = J_^?l _ 2C ^)9l ( 1 _ dK) ^ (3Q) 
6m* 2 dm q 3m* 3 y dm q J ' 



and 



24m* 2 dm q 1 2m* 3 y dm q ) 



We use <7tv = 45 MeV [25]. The derivatives of the scalar mean- field and the coupling 
constants with respect to the quark mass are calculated numerically - eg., the coupling 
constants are approximately given by quadratic or linear functions of the quark mass: 
for R = 0.8 fm, g 2 a = 269.6 - 2.150m + 0.01438m 2 , £ = 11.97 + 0.1616m , g 2 p = 
63.38 - 0.076m . 

Figure 9: The ratio of the quark condensate in medium to that in vacuum (m = 5 MeV 
and Rq = 0.8 fm). The solid and dotted curves are, respectively, for f p = 0.5 and 0. 



In Fig.|9| we show the ratio of the quark condensate in medium to that in vacuum. As 
one can see, the density dependence of the quark condensate can be well fitted by a linear 
function of the density : 

9kA „ ! _ (es\ , 32) 

Q(0) - 1 10.34JUJ' [ } 

where the reduction factors, (9 34), are for the proton fractions, f p = (° 5 ), respectively. 
The value of 0.36 agrees with the model-independent prediction of Cohen et al. for 
symmetric nuclear matter. Note that the ratio is not sensitive to the bag radius. If we 
take a heavier current quark mass (eg., m = 10 MeV), the reduction in the ratio becomes 
smaller (about 1 % at pg — 2p ) than the present case. 
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Taking up the terms of O(ps) in Eq.(^4j) to see the behavior at low density, the ratio 
is given by 

Q(Pb) 1 3ojv ( m a \ 

Q(0) C N (Q)m 2 J 2 \g a J 

Using Eq.(Trj), g 2 a ~ 260 and C N {0) ~ 0.37, Eq.® can reproduce Eq.fl32|). On the 
other hand, since the ratio of the hadron mass in medium to that in free space is roughly 
described by Eq. , the ratio of the quark condensates could be linked to the ratio of 
the hadron masses at low density: 



hadron ^ 



0.124 x n (l-^4). 



Mhadron \ Q(0) J 

This relation suggests that the ratio of the effective hadron mass to the free one does not 

1 /3 

always scale as ( ^§P ) II- At small density we could expect Eq . (|22l) , fl23T) and 

1/4 



(<\ „ (Q(pb)X 

\mj ~ ^ Q(0) J 

1 /3 

rather than the naive scaling. 



(35) 



V CONCLUSION 

We have applied the quark-meson coupling (QMC) model to investigate the density de- 
pendence of the properties of hadrons and the quark condensates in dense nuclear matter 
(up to four times the normal nuclear density). We have calculated not only the variations 
in the nucleon properties in medium but also those in the masses of the vector (u, p) 
mesons and the hyperons (A, S, S). As several authors have suggested 0, || || ]8[], the 
hadron mass is reduced due to the change of the scalar mean-field in medium. In the 
present model the hadron mass can be related to the number of non-strange quarks and 
the strength of the scalar mean-field. The hadron masses are simply connected to one an- 
other, and the relationship among them is given by Eq . (|22|) , which is effective over a wide 
range of the nuclear density. Furthermore, the ratio of the quark condensate in medium 
to that in vacuum can be related to the vector-meson mass and the nucleon properties - 
eg., g* A and p* N - in medium. 
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Finally, we would like to give a few caveats concerning the present calculation. The 
basic idea of the model is that the mesons are locally coupled to the quarks. This is 
certainly common for pions in models like the chiral or cloudy bag f26f , but may be less 
justified for heavier mesons like the vector mesons, which are obviously not collective 
states. The model also omits the effect of short-range correlations among the quarks. At 
very high density these would be expected to dominate. Furthermore, the pionic cloud 
of the hadron [^6| should be considered explicitly in any truly quantitative study of the 
hadron properties in medium. As we mentioned at the end of Sec. II, the contribution 
of meson exchange currents in medium is also important and should be considered in a 
consistent manner. Finally, we note that subtleties such as scalar- vector mixing in medium 
and the splitting between longitudinal and transverse masses of the vector mesons have 
been ignored in the present mean-field study. Although the former appears to be quite 
small in QHD the latter will certainly be important in any attempt to actually measure 
the mass shift. 

We would like to thank A. G. Williams for helpful comments on the manuscript. This 
work was supported by the Australian Research Council. 
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